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We have investigated properties of a-decay chains of recently produced superheavy elements 
Z = 115 and Z = 113 using the new Lagrangian model NL-SVl with inclusion of the vector self- 
coupling of LO meson in the framework of the relativistic mean-field theory. It is shown that the 
experimentally observed alpha-decay energies and half-lives are reproduced well by this Lagrangian 
model. Further calculations for the heavier elements with Z=117-125 show that these nuclei are 
superdeformed with a prolate shape in the ground state. A superdeformed shell-closure ai Z = 118 
lends an additional binding and an extra stability to nuclei in this region. Consequently, it is 
predicted that the corresponding Qa values provide a-decay half-lives for heavier superheavy nuclei 
within the experimentally feasible conditions. The results are compared with those of macroscopic- 
microscopic approaches. A perspective of the difference in shell effects amongst various approaches 
is presented and its consequences on superheavy nuclei are discussed. 

PACS numbers: 21.60.-n, 21.30.Fe, 27.90.-|-b, 23.60.-|-e 
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I. INTRODUCTION 

The large progress in superheavy element research has 
delivered challenging experimental data. This provides 
nuclear structure theory an opportunity to interpret the 
new experimental results and to improvise on predictions 
for the location of the superheavy nuclei on the basis of 
the new experimental data. At present there exist two 
sets of data due to two methods employed for superheavy- 
element synthesis 

The heaviest elements known today, the elements with 
Z = 110, 111, and 112, have been produced by cold fu- 
sion of heavy ions. Targets of ^"^Pb or ^°^Bi were irra- 
diated with the appropriate projectiles ^^Ni or ''"Zn to 
produce the isotopes ^"^^110, ^'^^111, and ^77112 0, 
These have been identified by long a decay chains lead- 
ing to known isotopes. These data are on a safer ground, 
for they have been verified in different laboratories^ 1^. 
Recently, even the heaviest species in this series, '^'^'^112, 
has been confirmed independently 0, ■ The important 
new result of these experiments is the discovery of a re- 
gion of deformed, shell stabilized nuclei 0, centered at 
Z = 108 (Hassium) and N ~ 162. Theoretical calcula- 
tions have argued that these nuclei are stabilized due to 
a hexadecapole deformation in the ground state |l.l^.ll0|. 

The second approach to synthesize heavy elements 
leads to elements with Z = 114, 115, and 116. Beams of 
^^Ca on actinide targets e.g. ^'*'*Pu and ^''^■^''^Cm were 
used. The isotopes of 286-289114^ ^^^^ 290-293115 ^^^^^ 

been identified in various experiments pH ITsL H^ . 
Recently, the isotopes 287,288115 i^^ye been observed in 
irradiations of ^^■^Am target with '*^Ca beam PI. These 
results are exciting and of special attraction for theory 
as the new region is located already close to the expected 
superheavy nuclei. The experimental problem is that the 



observed nuclides decay over long a-decay chains ending 
in spontaneous fission. These form an island of nuclei 
in itself and cannot be connected to the known region 
of isotopes. This makes an unambiguous identification 
with the presently used parent-daughter method impos- 
sible. To cope with this problem, a number of consistency 
checks have been made (14| . In addition, first promising 
experiments to verify the chemical nature of these ele- 
ments are underway p^ . 

Experimental investigations of superheavy elements 
encounter a general problem that only a few atomic nuclei 
of each species are produced. Presently, for the heaviest 
species only basic properties can be directly extracted 
from experiment. These are Qa-values and lifetimes. 
Most of the known nuclides have odd proton or neutron 
numbers and are not likely to decay by ground state tran- 
sitions. Moreover, superheavy nuclei have a complicated 
level structure with high spin and low spin states be- 
ing close together. They are likely to form isomers 
This makes a comparison with theoretical predictions a 
little difficult. In those cases where the chains end in fis- 
sion, Qa-values and lifetimes are the only experimental 
information that can be extracted for nuclides along the 
chain. However, ground state masses of nuclei cannot be 
obtained. 

Properties and structure of superheavy nuclei have 
been investigated extensively using various approaches. 
The approaches consist of microscopic nature such as 
non-relativistic density-dependent Skyrme Hartree-Fock 
(SHF) theory and the Relativistic Mean-Field (RMF) 
theory or macroscopic-microscopic type. In the latter 
category, total binding energy of nuclei is obtained as 
a sum of a smooth energy based upon liquid drop type 
formula on which shell correction is imposed using the 
method of Strutinsky pTlll^ . The most notable effort in 
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this direction has been the Finite-Range Droplet Model 
(FRDM) The shell correction energies were cal- 

culated in the FRDM in order to identify major magic 
numbers in the region of superheavy nuclei. The FRDM 
predicts a major proton magic number at Z = 114 be- 
yond the well known magic number Z = 82. Experi- 
mental data, however, give little support to this magic 
number. Calculations using the macroscopic-microscopic 
approach have also been performed by another group 
[20I I21I |23 |. where Yukawa-plus-exponential model 
has been used for the macroscopic component and the 
Strutinsky shell correction is used for the microscopic 
component. Calculations of properties of superheavy 
nuclei have been rather successful in the macroscopic- 
microscopic approaches. 

Extensive studies of superheavy nuclei have also been 
performed within self-consistent mean-field models of 
SHF type [23| or c omp arative studies of SHF models and 
RMF models [H 113, El 113 have been done. One 
of the first studies of the properties of superheavy nu- 
clei within the RMF theory was performed earlier in ref. 
[30I I using the force NL-SH 31]. It was shown that neu- 
tron number = 184 appears to be magic for lighter 
superheavy nuclei, whereas it diminishes for heavier su- 
perheavies for Z > 114. In this work shell correction 
energies were also calculated for heaviest deformed su- 
perheavy nuclei. For the proton number, there was an 
indication of a deformed shell closure at Z = 114. This 
was consistent with the predictions of the FRDM for a 
strong shell gap at Z = 114 in the deformed region. 

In order to identify magic numbers in the region of 
superheavy elements theoretically, extensive shell correc- 
tion calculations have been performed by Kruppa et al. 

for a set of the Skyrme and RMF forces. It is shown 
that Skyrme models predict the strongest shell effects at 

= 184 and Z = 124, 126, but not at Z = 114. On the 
other hand, several RMF forces considered in this work 
[s^l do not show a shell gap N — 184. This contrast- 
ing difference in the behavior of magic numbers between 
the Skyrme forces and RMF forces can be attributed to 
the difference in the shell structure in the RMF and SHF 
approaches. This derives partly from the difference in 
the spin-orbit interaction and its isospin dependence. It 
was shown in ref. [33j that isospin dependence in the 
Skyrme Hartree-Fock theory is different than that in the 
RMF theory. This difference in the isospin dependence 
of the spin-orbit potential was shown to be responsible 
for a correct description of anomalous isotope shifts in 
Pb isotopes '33|. A modification in the isospin depen- 
dence of the spin-orbit potential ,33> .34] in the Skyrme 
Ansatz allowed a proper description of the isotope shifts 
in Pb nuclei, though refs. [s^] and provide for dif- 
ferent prescriptions for the dependence. It is, therefore, 
expected that Skyrme forces without such a modification 
would yield results different than those from the RMF 
forces. 

A significant difference between the RMF theory and 
SHF theory has emerged in a systematic study of fission 



barriers of actinide nuclei. A comparative study of fission 
barriers has been performed recently in SHF models and 
RMF theory j^i]. It has been shown that the RMF forces 
considered in this work predict lower fission barriers than 
most of the SHF models. It is, however, not clear whether 
this discrepancy holds for all the available RMF forces. 
This difference between RMF and SHF forces, if it holds 
in general, needs to be understood. It is surmised that 
this difference may lie in the difference in shell structure 
between the two theories. 

We have studied the recent a-decay chains of Z = 115 
[THI ] employing the RMF theory in the present paper. As 
it is well established, the Relativistic Mean-Field (RMF) 
theory has proved to be successful in providing a 
framework for description of various facets of nuclear 
properties [11113 EE l22l|4ii|4ll. The relativistic Lorentz 
covariance of the theory allows an intrinsic spin-orbit in- 
teraction based upon exchange of a- and w-mesons. This 
has been shown to be advantageous for properties such 
as anomalous isotope shifts in Pb nuclei [42] ■ The RMF 
theory has achieved significant success in respect of nuclei 
near the stability line as well as for nuclei far away from 
the stability line |31., jiS^ ^] ■ The nonlinear scalar self- 
coupling of cr-meson has been the most successful model 
in the RMF theory used so far. However, a closer anal- 
ysis of Lagrangians with the non- linear scalar coupling 
of cr-meson has shown (4^ that shell effects with these 
model Lagrangians at the stability line are stronger than 
the experimental data. This effect is extended to regions 
far beyond the stability line. It has been shown that the 
nature of the shell effects along the stability line does in 
turn influence the shell effects in the extreme and un- 
known regions jil ]. 

In the present paper, we have calculated properties of 
Qf-decay chains of the newly discovered superheavy nuclei 
^^^115 and ^^^115 jl3] using the Lagrangian model of the 
vector self-coupling of uj meson. In Section II, we present 
a perspective of the shell effects with a view to show a 
relationship of the shell effects in the r-process region to 
those in the superheavy region. We allude briefly to some 
of the distinctive features of the RMF theory in Section 
II. We will also refer to some of these features in our 
discussion. A brief formalism of the RMF theory is pre- 
sented in Section HI, and details of RMF calculations are 
provided in Section IV. Section V discusses results of the 
RMF calculations on Qa values, a-decay half-lives and 
deformation properties. A comparison of the results on 
Qa values and a-decay half-lives is made with some ear- 
lier calculations and also with predictions made by some 
macroscopic-microscopic models. We also predict prop- 
erties of heavier superheavy nuclei which can be a-decay 
precursors of the nuclei ^^*115 and ^^^115. Deformed 
single-particle spectra are presented in order to visualize 
existence of possible islands of stability in the deformed 
region. The last section summarizes the results in view 
of experimental feasibility of synthesis of heavier super- 
heavy nuclei. 
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II. IMPORTANCE OF THE SHELL EFFECTS 

The spin-orbit interaction and consequently how the 
shell effects behave in the extreme regions plays a sig- 
nificant role in carving out shell gaps in superheavy nu- 
clei and in nuclei near r-process path. In this respect, 
the RMF theory has an inherent advantage in that the 
spin-orbit interaction arises naturally as a result of the 
Lorentz-Dirac structure of nucleons. Consequently, the 
RMF theory has shown an immense potential in being 
able to describe properties of nuclei along the stability 
line and also for a large number of nuclei beyond the 
stability line. 

Recently, it was shown ,43j that Lagrangian models 
with the non-linear scalar coupling of a meson only over- 
estimate the shell effects at the stability line. In order to 
remedy the problem of strong shell effects, the nonlinear 
vector self-coupling of uj meson was introduced in Ref. 
[43l| . It was shown ^3 that the ensuing Lagrangian pa- 
rameter set NL-SVl is able to reprodduce the shell effects 
in Ni and Sn isotopes near the stability line. It was also 
shown j4^] that the force NL-SVl based upon the vector 
self-coupling of w-meson is also successful in describing 
the available data on the shell effects across the waiting- 
point nucleus *'^Zn. This exemplifies the importance of 
an appropriate description of the shell effects along the 
line of stability with a view to be able to extrapolate 
the shell effects in the extreme regions such as near the 
r-process path. Due to this reason, we consider that a 
description of superheavy nuclei in the extreme regions 
can be put on a footing similar to that of nuclei which 
are close to the r-process path. 



TABLE I: The Lagrangian parameters of the force NL-SVl 
used in the calculations. All the masses are in MeV. While g2 
is in fm~^, all the other coupling constants are dimensionless. 
Here g4 represents the quartic coupling of ui meson. 



M = 939.0 




=10.1248 


92 = 


-- -9.2406 


m<,=510.0349 




=12.7266 


S3 = 


-15.388 


m„=783.0 


9p 


= 4.4920 


94 = 


41.0102 


mp=763.0 











scalar potential of the form: 

U{a) = \mW + + \9^o\ (2) 

where 172 and 33 represent the parameters of the nonlinear 
coupling of cr meson to nucleon. This model Lagrangian 
has been very successful in describing properties of nu- 
clei at the stability line and also for nuclei away from the 
stability line. The coupling constant for the non-linear lo- 
term is denoted by 54 in the Lagrangian of Eq. (1). The 
vector self-coupling of the w-meson was first proposed 
in reL pj||, where properties of nuclear matter associ- 
ated to this potential were discussed. It is interesting to 
note that introduction of the non-linear coupling of lo- 
meson softens the equation of state (EOS) of the nuclear 
matter significantly. This has the consequence that the 
maximum neutron star mass with such an EOS appears 
within the bounds of empirically observed values. 



IV. DETAILS OF THE CALCULATIONS 



III. RELATIVISTIC MEAN-FIELD THEORY 

The RMF approach 36] is based upon the Lagrangian 
density which consists of fields due to the various mesons 
interacting with nucleons. The mesons include the 
isoscalar scalar cr-meson, the isovector vector w-meson 
and the isovector vector p-meson. The Lagrangian den- 
sity is given by: 

C = ^{i)- g^i^ - gppT - ie(l - T3)4 - g^a - ?/) 

^\d^,ad^u - U{a) - ^n^M'^'^ + ^mluj^ajf^ (1) 

The bold-faced letters indicate the vector quantities. 
Here M, m^., and uip denote the nucleon-, the cr-, 
the to- and the p-meson masses respectively, whereas go-, 
guj, gp and e^/47r — 1/137 are the corresponding couphng 
constants for the mesons and the photon, respectively. 

The nonlinear scalar interaction has been used exten- 
sively to describe the ground-state properties of finite 
nuclei. Herein, the a meson is assumed to move in a 



The RMF calculations in this work have been per- 
formed in an axially symmetric deformed basis. The 
method of expansion |3q of wavefunctions into harmonic- 
oscillator basis has been used to solve the Dirac and 
Klein-Gordan equations. Both the fermionic and bosonic 
fields have been expanded in a harmonic oscillator ba- 
sis of 20 shells. The pairing has been included within 
the BCS scheme employing constant pairing gaps. The 
prescription of reL ^j] which provides for a best fit of 
pairing gaps for neutrons and protons over a large range 
of nuclei has been used. Accordingly, pairing gaps are 
taken as 

A„(,) = 4.87V-i/3(Z-i/3) (3) 

where N and Z represent the neutron and proton number 
of a nucleus. The Lagrangian parameter set NL-SVl 
has been used for all the calculations. The parameters 
of the force NL-SVl are presented in Table I. In this 
work, we have not included the blocking of pairing for an 
odd particle. Extensive calculations with the blocking of 
pairing would be presented elsewhere. 

We have carried out constrained calculations with a 
quadratic constraint. The potential energy landscapes 
of nuclei over a large range of quadrupole deformation 
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have been mapped out and various possible minima in 
the total energy have been identified. This also serves 
to confirm the total binding energy of ground state as 
obtained in individual RMF+BCS minimizations. 



V. RESULTS AND DISCUSSION 

A. Qa values and half-lives 

The decay chain of the superheavy nucleus ^^^115 
(Z=115) was observed |0 for 3 events in the 3n evap- 
oration channel in the reaction ^''^Am + ^®Ca with the 
projectile energy E = 248 MeV. The energy of alpha 
particles emitted as a result of consecutive a-decay of 
various product nuclei is within about 250 keV in the 
three events. The maximum dispersion in the alpha en- 
ergy has been observed for ^^^113 amounting to about 
0.5 MeV. The best two values agreeing with each other 
have been selected as the acceptable value for Qa in ref. 
Il5|. 

Results of the RMF calculations with the Lagrangian 
set NL-SVl for the newly discovered chains of superheavy 
nuclei are shown in Figs. and [3 Figure ^ shows the 
difference AQa of the calculated Qa values from the ex- 
perimentally obtained ones for the chain with the parent 
nucleus ^^^115. For the nucleus ^^^115, calculations with 
NL-SVl show a deviation from the experimental value by 
about 0.8 MeV. For the nuclei ^^^113 and ^^^111, NL-SVl 
shows a very good agreement with the data. On going to 
lighter products of the chain, i.e. for ^^^Mt (Z=109) and 
^^^Bh (Z=107), our results agree with the data within 
about half an MeV. A comparison is also made with the 
relativistic mean-field calculations of ref. using the 
parameter set TMA HP]. Results of TMA with the use of 
constant pairing (open square) and a density-dependent 
(5-pairing force (diamonds) are shown. The results shown 
with the (^-pairing also include the blocking effect. The 
Qa obtained with the constant pairing for the nucleus 
^^^115 is about the same as that with NL-SVl. However, 
for other nuclei such as 284ii3 and 272Bh (Z=107) the 
disagreements with the data amounts to about 1 MeV. 
Whilst for ^^''113, TMA with constant pairing overesti- 
mates the experimental value by about 1 MeV, it un- 
derestimates the experimental Qa for ^^^Bh (Z=107) by 
about the same amount. In comparison, the use of the 
i5-pairing force with the blocking has shown some im- 
provements in the results with TMA. It is seen from the 
figure that our present calculations with NL-SVl are in 
good agreement with the Qa values of the decay chain of 
^^^115 with the exception for the parent nucleus. 

In ref. a single decay chain of ^^^115 was observed 
on increasing the projectile energy by 5 MeV to E = 253 
MeV. It is assumed that 4n evaporation channel in the 
excitation function becomes dominant. It may, however, 
be said that any other event for this decay chain has not 
been observed to confirm the a-decay energies. In view 
of this, these values may be treated with caution. The 
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FIG. 1: The difference AQa — Qa (theory) - Qa (expt.) be- 
tween the theoretical and the experimentally observed Qa val- 
ues for the a-decay chain with the parent nucleus '^^^IIS. The 
RMF results with NL-SVl (full circles) are compared to those 
with the RMF calculations using the force TMA with constant 
pairing gap (open squares) and using a density-dependent 5- 
pairing force including blocking effect (open diamonds). 
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109 111 113 115 

Z 

FIG. 2: The same as for Fig. Q but for the Q-decay chain 
with the parent nucleus ^^""^115. 
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FIG. 3: The Qq values obtained with the force NL-SVl in 
the RMF calculations are compared with the experimental 
data from ref. for the a-decay chain of ^^^115 in the 

upper panel. The Qa values predicted by our model for the 
heavier parents of this chain are also shown. In the lower 
panel the ground-state quadrupole deformation [32 are shown 
by open circles. The solid circle for the nucleus ^^^113 exhibits 
a second minimum at the highly deformed (superdeformed) 
configuration that is in coexistence with the low deformed 
prolate state. 



a-decay energies for the decay chain of ^^'115 as reported 
in the ref. ^ are 10.59 MeV, 10.12 MeV, 10.37 MeV 
and 10.33 MeV, respectively. The two latter values show 
an increase over that of the second value (10.12 MeV) in 
the decay chain. This behaviour is different from that of 
the ^^'^115 decay chain, where the a-decay energies show 
a decreasing trend with the decrease in the Z value of the 
daughter nuclei. Assuming that the shell structure has 
not changed from the ^^^115 chain to the ^^^115 chain, 
there could easily be an uncertainty of about 0.5 MeV in 
the a-decay energies of ^^^111 and ^^^109. 

The results of our calculations for the decay chain of 
28^115 are shown in Fig. 2. The values with NL-SVl 
agree well with the experimental values for the nuclei 
^^''115 and ^^^113. For the nucleus ^''^111, our calcula- 
tions show a Qa of 10.0 MeV compared to 10.52 MeV 
as observed in ref. J^] . The difference is accentuated to 
about 1 MeV for the nucleus ^"^^109. It is noted that the 
calculated values with NL-SVl show a decreasing trend 
with the Z value. An uncertainty in the a-decay energies 
for this decay chain as commented above may partly ac- 
count for this discrepancy. Generally, there is an overall 



FIG. 4: The same as for Fig. 3 but for the chain consisting 
of the nucleus ^**^115. 



agreement of the NL-SVl results within 0.5 MeV for the 
Qa values save that for ^'^^Mt (Z = 109). 

Calculations with TMA from ref. 4S! are also com- 
pared in Fig. 2. The calculations both with constant 
pairing and with (5-pairing with blocking show a trend 
similar to that of NL-SVl vis-a-vis the experimental val- 
ues. In the calculations of ref. Ji^l disagreements are 
more pronounced for ^ssna and ^^SMt {Z = 109). 

B. Predictions for heavier superheavy nuclei 

We have calculated ground-state properties of the de- 
cay chains of ref. . We have also extended our calcu- 
lations in order to be able to predict properties of odd-Z 
nuclei heavier than ^^^115 and ^^"^115, whereby ^^^115 
and ^®^115 would constitute daughter nuclei of a sequen- 
tial a-decay chain. The Qa value and quadrupole de- 
formation (32 corresponding to the lowest energy state 
for the chain comprising the ^^^115 and its daughters are 
shown in Fig. 3. The NL-SVl results show general agree- 
ment with the experimental Qa values !lj^ (see Table II) 
as discussed above and shows an increasing trend with 
an increase in the Z value. The NL-SVl value, on the 
other hand, shows a slight decrease at Z = 115. This is 
indicative of a possible deformed shell closure at Z = 114 
with NL-SVl. For all the nuclei above Z = 115, the Qa 
value shows a systematic increase with a slight kink at 
Z = 119. 

The corresponding (32 values are shown in the lower 
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TABLE II: The Qa values, a-decay half- lives Ta, and quadrupole deformation /32 of nuclei in the decay chain consisting of the 
nuclide 2**115 as obtained with RMF calculations using the Lagrangian set NL-SVl. The values and estimated half-lives 
from the recent experimental results [l5|| are also shown for comparison. 



Nucleus 


Qc (MeV) 




P2 


(exp.) (Me 


308^25 


13.41 


16 /is 


0.62 




304^23 


12.64 


0.20 ms 


0.61 




300 -1^21 


11.74 


6.5 ms 


0.56 




296119 


10.81 


0.38 s 


0.54 




292117 


10.35 


1.63 s 


0.52 




2**115 


9.79 


15.3 s 


0.50 


10.61 ±0.06 


2*^13 


10.27 


0.15 s 


0.17 


10.15 ±0.06 


2*0lll 


9.75 


0.97 s 


0.18 


9.87 ± 0.06 


2^6 Mt 


9.40 


2.21 s 


0.20 


9.85 ±0.06 


2^2Bh 


8.60 


150 s 


0.26 


9.15 ±0.06 



Tc (exp.) 



i°«ms 

+0.58 
0.17* 

3.6t*-3 



87 
0.48 



0.72 



1.3 ° 

+0.87 
-0.25 
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o 
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■a 
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Z 
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_i_ 
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A 



FIG. 5: A comparison of the Qa values obtained with NL- 
SVl is made with the resuhs of ref. (2^ and FRDM [HI 
for the chain comprising 2**115. The experimental values 
obtained from the recent experiment are also shown for 
comparison. 



panel. The quadrupole deformation stays in the vicinity 
of ^2 ^ 0.2 in going from Z = 107 to Z = 113. How- 
ever, the nucleus ^^"^IIS marks a transition point for the 
chain, in that it exhibits a shape coexistence of a low 
prolate deformation [Ji^ — 0.17) with a high prolate de- 
formation (/32 — 0.50) in the ground state. Here, the 
state with the low deformation is the lowest one. All the 
nuclei including and above Z = 115 exhibit a large pro- 
late deformation {(32 ^ 0.6) in the ground state. This is 
akin to a superdeformation. As the single-particle scheme 



FIG. 6: The same as for Fig. 5 but for the chain consisting 
of the nuclide 2*^115. 



would show, this is associated with a strong shell gap at 
Z = 118 in the superdeformed region. This allows an 
increase in the binding energies of nuclei and accordingly 
the Qa values are comparatively suppressed. 

The Qa and (32 values for the ground state of nuclei 
in the chain associated with ^^^115 are shown in Fig. 4. 
As discussed above, the Qa value for the nucleus ^^^Mt 
shows a deviation of about 1 MeV from the observed 
value. For the others, there is a good agreement with the 
experimental value (see Table HI). The NL-SVl values 
show a systematic increase in Qa in going from Z = 109 
to Z = 115. A kink in Qa at Z = 115 indicates a presence 
of deformed shell closure at Z = 114. 
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TABLE III: The same as for Table II but for the decay chain consisting of the nuclide 115. 



Nucleus 


(MeV) 




P2 


Qc (exp.) (MeV) 


r„ (exp.) 


^•"125 


13.51 


10 us 


0.61 






303;^23 


12.86 


68 fis 


0.60 






299121 


11.91 


2.6 ms 


0.55 






295119 


10.94 


0.17 s 


0.54 






291117 


10.26 


2.9 s 


0.52 






287^-^5 


10.57 


100 ms 


0.51 


10.74 ±0.06 


32+if ms 


283113 


10.35 


90 ms 


0.17 


10.26 ±0.06 


WOttf ms 


2^^111 


10.00 


200 ms 


0.19 


10.52 ±0.06 


170+80° ms 


2^5 Mt 


9.43 


1.8 s 


0.21 


10.48 ± 0.06 


9.7t« ms 



As far as deformation of nuclei in this chain is con- 
cerned, all nuclei with low Z values including ^^'^113 ex- 
hibit a low value of /32 in the vicinity of 0.20. Here, the 
nucleus ^^^113 with Z = 113 exhibits the phenomenon 
of shape coexistence of a low prolate deformation {P2 — 
0.17) with a high prolate deformation (/?2 — 0.50) in the 
ground state. The state with with the low (32 value is the 
lowest minimum in the binding energy. This behavior is 
similar to that exhibited by the corresponding nucleus 
^*'*113 in the other chain. Likewise, the nucleus ^^'^113 
marks a transition point in the shapes of nuclei in this 
chain. All the nuclei including and above Z = lib show 
a superdeformed ground state. All the features shown in 
Fig. 4 for this chain exhibit a similarity to those shown 
in Fig. 3. Therefore, properties of nuclei for both the 
chains point to almost the same structural behavior. 

RMF calculations with another Lagrangian parameter 
set NL-SV2 [i^, with the vector self-coupling of iv- 
meson for both the chains of Figs. 3 and 4 show that 
nuclei ^*'*113 and ^^^113 exhibit a shape coexistence of 
low prolate deformation with a high prolate deformation 
in the ground state. The magnitude of deformations with 
NL-SV2 is seen to be similar to that obtained with NL- 
SVl as shown in Figs. 3 and 4. With the Lagrangian 
set NL-SV2, superheavy nuclei heavier than ^^^113 and 
^^■^113 in these chains are also observed to be superde- 
formed in the ground state. Accordingly, the nuclei 
^^'*113 and ^^^113 denote a transition point from a low 
deformation to a superdeformation also with NL-SV2. 
Thus, the results on the shape-coexistence and the pres- 
ence of a second superdeformed minimum for nuclei with 
Z = 113 with NL-SVl are also exhibited by NL-SV2. A 
comprehensive comparison of results in the superheavy 
region wtih various Lagragian models will be presented 
elsewhere [s^. 

In a recent investigation of deformation properties 
of superheavy nuclei it has been noted that 

macroscopic-microscopic approach does not support a su- 
perdeformation in the superheavy region. This differ- 
ence in the outcome of the results of ref. |51|] can be 
understood due to strong shell effects prevalent in the 
macroscopic-microscopic approaches. The strong shell 
effects or equivalently strong shell gaps would disfavour 
a superdeformation. 



Our results for the recently discovered chain of ^**115 
and predictions for heavier superheavy nuclei are com- 
pared with some other models in Fig. 5. Detailed 
macroscopic-microscopic calculations for a large number 
of nuclei in the superheavy region of Z — 110 — 120 have 
recently been performed by Muntian et al. The 
results available on some nuclei in this chain are shown. 
Whilst this model shows a reasonably good agreement 
with the experimental values on ^^^115 and ^^^113, it 
overestimates the datum for ^*°111. It is interesting to 
compare the data also with those from the macroscopic- 
microscopic model FRDM js^. The FRDM values show 
an agreement with the data for two lighter superheavy 
nuclei and also for ^^^115 within about 0.5 MeV. How- 
ever, for the nucleus ^^'^113 the FRDM value underesti- 
mates the Qa value by ~ 1 MeV. For heavier a-decay 
precursors of ^^^115, i.e., nuclei with Z = 117 — 125, we 
can compare our results with NL-SVl with those from 
FRDM. 

As far as a comparison of the macroscopic-microscopic 
calculations of Muntian et al. with FRDM are con- 
cerned, the predictions only for ^^^117 can be compared. 
For this nucleus the two models agree. As further results 
on other odd nuclei from ref. are not available, it is 
difhcult to compare the two models in the context of the 
two chains studied in our work. However, as predictions 
from ref. are generally on the higher side in this re- 
gion, we surmise that the results of the two macroscopic- 
microscopic would follow each other for heavier super- 
heavies. In comparison, our results with NL-SVl show a 
smooth increase in the value of Qa with Z value. The 
NL-SVl results, on the other hand, predict Qa values 
which are smaller by about 1-2 MeV from those of the 
FRDM. This is due to an extra stability lent by the su- 
perdeformed ground state of all the heavier nuclei with 
Z > 115. This will be demonstrated by a large shell gap 
in the deformed single-particle spectrum at Z = 118 and 
Z = 120. 

Calculation of half-lives of a-decay from Qa values 
is not free of uncertainties. This is subject to Viola- 
Seaborg systematics fitted over a given region of nuclei. 
The mostly used systematics derive from a fit of param- 
eters done in ref. 53:]. We have calculated the a-decay 
half-lives of nuclei using the formula from Viola-Seaborg 
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systematics. Accordingly, the a-decay half-life can be 
written as 

log T^{s) = (aZ + b)Q-^^^ + (cZ + d) (4) 

where Z is proton number and Qa is alpha-decay energy 
of the parent nucleus in MeV. The parameters are a — 
1.66175, b = -8.5166, c = -0.20228 and d = -33.9069. 
These parameters are taken from ref. '5^, where these 
were readjusted in order to take into account new data. 
These parameters have been found to be successful over 
a broad range of nuclei. 

The Ta values calculated using Qa values from NL- 
SVl for nuclei of the chain consisting of ^^^115 are shown 
in Table II. The experimentally observed a-decay half- 
life ^3 of ^^^^115 and its daughters are also shown for 
comparison. It may be noted that experimental values 
of Ta for this chain are based upon 3 events observed in 
the experiment [l5l |. The corresponding half- life quoted 
in ref. TE'l has been obtained from the average of the 
life-times observed in the experiment. Consequently, one 
can notice large error bars in the experimental values. 
These values therefore signify the order of magnitude of 
a Ta value. 

Given the large uncertainties in the experimental val- 
ues, Ta values calculated with NL-SVl agree well with 
the experimental ones with the exception for ^^*115 and 
^''^Bh (see Table II). The theoretical half-life for these 
two nuclei are about an order of magnitude larger than 
the experimental values. This is due to the reason that 
with our model the Qa values are ~ 0.50 — 0.80 MeV 
smaller than the experimental Qa value. For other su- 
perheavy nuclei with Z = 109 — 113, there is a good 
agreement between the theoretical and experimental a- 
decay half-lives. This is due to a good agreement between 
the corresponding Qa values. 

In going to the heavier precursor nuclei in the same 
chain (Table II), there is an increasing tendency in the 
value of Qa with an increase in the Z value. However, Qa 
values with NL-SVl are systematically lower than those 
of FRDM. The lower Qa values with NL-SVl has the 
consequence that superheavy nuclei with high Z such as 
Z = 117, Z = 119 and Z = 121 would not have a very 
small a-decay half-life as compared to that from FRDM. 
The a-decay half-lives predicted from NL-SVl for the 
nuclei 292117, ^^^119, and ^00121 are - 1.6 s, 0.4 s and 
7 ms, respectively. This would put these nuclei within 
the range of experimental accessibility and feasibility. 

A comparison of various results for the chain compris- 
ing 287j^j^5 g^j-g shown in Fig. 6. In general, the results 
of Muntian et al. overestimate the new data by ^ 
0.5-1.0 MeV. On the other hand, the results of FRDM 
for this chain are very similar to those for the other chain 
as shown in Fig. 5. Though the FRDM results show a 
good agreement with a few data, kinks at Z = 111 and 
Z — 113 show an undulating character of the FRDM re- 
sults. Such kinks are not present in other results. This 
may be due to a strong spherical shell closure at Z = 114 
as predicted by the FRDM. 



The a-decay half-lives for the decay chain of 2®'115 
from NL-SVl are compared with the experimental esti- 
mates from ref. in Table III. It may be noted that the 
experimental value is based upon life-time estimates from 
a single decay event of the superheavy nucleus ^^^115. 
There is an overall good agreement between NL-SVl val- 
ues and the experimental ones. Only for ^^^Mt {Z = 109) 
does NL-SVl overestimate the half-life by about 2 orders 
of magnitude. The experimental half-life for the nucleus 
^''^Mt has been estimated to be 9 . This value 
is much lower than the experimental estimates for the 
heavier superheavies such as ^^^111 and ^^^113 (see Table 
III). In comparison, the experimental Qa value for ^^^Mt 
is comparable to that for ^^^111. According to Viola- 
Seaborg formula, this experimental Qa for ^^^Mt would 
indicate Ta value much larger than the quoted one |l5j . 
In view of this, there seems to be a little discrepancy be- 
tween the experimentally deduced Qa and Ta value for 
the nucleus ^^SMt. This may be due to ambiguities in 
identification in a single event. 

The generally higher Qa values from macroscopic- 
microscopic calculations 22] vis-a-vis the experimental 
data [l^ by ^ 0.5 — 1.0 MeV would mean a decrease 
in a-decay half-life by about 1-3 orders of magnitude in 
the macroscopic- microscopic calculations However, 
it is worth mentioning that the results of ref. '22^ on a- 
decay energies describe well the data on and variation in 
the Qa values with neutron numbers for isotopic chains 
below Z = 109, as shown clearly in Fig. 2 of ref. [T5|. 

For the heavier superheavies of the chain comprising 
of the nucleus 2^'' 115 (Fig. 6) the FRDM values show a 
strong increase in the Qa value. This would make the 
half-life of some of these nuclei in the range of nanosec- 
ond or even smaller. In comparison, our predictions with 
NL-SVl show only a modest increase in the Qa values 
for nuclei heavier than Z = 115 (see Table HI). This 
behaviour is very similar to that exhibited by the other 
chain in Fig. 5. The a-decay half- life using the Viola- 
Seaborg systematics for nuclei ^^^117, '^^^119, and ^°^121 
in our model are 2.9 s, 0.17 s and 2.6 ms, respectively. 
These values are well within the reach of experimental 
feasibility. Thus, the NL-SVl results shown in Fig. 5 
and Fig. 6 puts several heavy elements within the range 
of feasibility of realization. 



C. Potential energy landscapes 

We have performed constrained calculations with a 
quadratic constraint in order to map out potential en- 
ergy landscape. The total energy of nuclei as a function 
of quadrupole deformation /?2 is shown in Figs. 7 and 8 
for superheavy nuclei in the two isotopic chains studied 
in this work. All the nuclei of chain comprising of both 
288115 (Fig. 7) and "^^^115 (Fig. 8) exhibit a prolate de- 
formation for all the nuclei from Z — 111 to Z = 125. 
A configuration in oblate space seems to be forbidden 
as suggested by the total potential energy curve for the 
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FIG. 7: Potential energy landseapcs of nuclei in the decay chain comprising ^^^115 for (a) lighter superheavies below Z < 117. 
The shape coexistence of a deformed and a superdeformed prolate shape is shown for ^®^113 in the inset. (b)Potential energy 
landscapes of heavy superheavies with Z > 119. All the nuclei show a superdeformed ground state as depicted in the inset for 
296119. 
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FIG. 9: The proton Nilsson single-particle levels obtained FIG. 10: The neutron Nilsson single-particle levels obtained 
from the RMF calculations with NL-SVl for low-deformed from the RMF calculations with NL-SVl for low-deformed 
prolate-shaped nuclei of the decay chain of ^**115 are shown. prolate-shaped nuclei of the decay chain of ^^*115 are shown. 



D. Single- particle levels 



superheavy nuclei shown in Figs. 7 and 8. 

Consistent with previous description, the nuclei ^^°111 
and ^^^111 are prolate deformed with a low value of 
quadrupole deformation. On the other hand, the nu- 
clei 284113 (Fig. 7(a))and ^^^113 (Fig. 8(a)) exhibit two 
minima, as discussed earlier. The two minima, which co- 
exist in energy are depicted in the insets of Fig. 7(a) and 
Fig. 8(a). The second minimum at P2 ~ 0.50 lies within a 
few hundred keV of the first minimum. For all the heav- 
ier elements with Z > 113, the ground state minimum 
is obtained in the region of a high prolate deformation 
as seen by a single well in the curves in part (a) and (b) 
of Figs. 7 and 8. The high prolate deformation of su- 
perheavies continues in the high Z region as exemplified 
in the insets of Fig. 7(b) and Fig. 8(b) for the element 
Z = 119. It is interesting to note that the potential en- 
ergy landscapes for Fig. 7 are very similar to those of 
Fig. 8. 



The Nilsson single-particle levels obtained with NL- 
SVl for the nuclei of chain ^^8115 ^re shown in Figs. 9 
and 10. The proton single-particle levels for nuclei in 
this chain are presented in Fig. 9. The corresponding 
ground-state (32 values for each nucleus are indicated in 
the top panel of the figure. It shows a decrease from 0.26 
to 0.17 in going from ^^^Dh {Z = 109) to 284ii3. This 
decrease is accompanied by an arrival of a deformed shell 
gap at Z = 114 as seen in the figure. There are also 
minor deformed shell gaps at Z = 106 and Z = 108. As 
the deformation is changing only slightly and the Fermi 
energy is showing a marginal change in going from the 
left to the right, the level crossing in Fig. 9 is not such 
a dominant feature. However, shifting of the deformed 
levels partly due to the change in deformation and partly 
due to the change in the chemical potential, deformed 
shell gaps are being created at Z — 108 for the lower Z 
nuclei and at Z = 114 for the higher Z nuclei. There is, 
however, no spherical Z = 1\A closed shell as suggested 
in models such as FRDM [ll|. 
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FIG. 11: The proton Nilsson single-particle levels obtained 
from RMF calculations with NL-SVl for superdeformed pro- 
late nuclei, which are a-decay precursors of the nucleus 
28«115. 



FIG. 12: The neutron Nilsson single-particle levels obtained 
from RMF calculations with NL-SVl for superdeformed pro- 
late nuclei, which are a-decay precursors of the nucleus 
28«115. 



The neutron single-particle levels in Fig. 10 show a 
shell gap at iV = 154 and N = 164. The latter is in- 
evitably associated with the nuclei ^^'^Db and ^^^Bh. For 
the heavier nuclei with A — 280 and A = 284 the de- 
formed shell gap evolves at N = 174. This has been 
observed to the case in several other models. It can be 
seen that a frequent crossing of levels creates a bunch- 
ing and thus deformed shell gaps at = 154, A'' = 164, 

= 174 and TV = 180 are produced. The well-known 
shell gap at = 184 is not seen in this figure, though it 
is not relevant to mass numbers shown in this figure. 

The proton Nilsson orbitals for the precursors of the 
nucleus ^*^115 are shown in Fig. 11. For these nuclei the 
ground state quadrupole deformation (32 increases from 
0.50 for 288115 to 0.62 for ^°^125. Because these nuclei 
are superdeformed, there is a significant level crossing 
in this diagram. A deformed shell gap ai Z = 112 is 
seen. However, this shell gap is much less relevant to the 
higher Z values of nuclei considered in this fig ure. For 
the experimentally observed nucleus ^^^115 [13, there is 
no perceptible observation of a deformed shell gap at Z = 
114. This gap which was seen in the lighter counterparts 



of the chain is washed out here. In contrast, a strong 
shell gap is observed at Z = 118. This shell gap in the 
highly deformed region provides an extra stability to the 
nuclei concerned. This has already been seen in terms 
of lower Qa values for nuclei with higher Z as shown in 
Fig. 5. A possible deformed shell closure at Z = 120 is 
also observed, but for very heavy superheavies. Here the 
nuclide ^°*125 runs close to the proton drip line. 

The corresponding neutron Nilsson orbitals for the su- 
perdeformed superheavy nuclei which are heavier parents 
of the nucleus ^^^115 are shown in Fig. 12. The neutron 
number for these nuclei spans the range A^ — 173 — 183. 
For the lighter counterparts of the chain one can see a 
deformed shell gap at A^ = 174, that diminishes in going 
to the heavier nuclei. There is another deformed shell 
gap at N = 186. However, we do not see a shell gap at 
A^ = 184 in the deformed region. Moreover, nuclei in the 
vicinity of this neutron number are assuming superde- 
formed prolate shapes. This discounts the possibility of 
a deformed shell gap at A^ = 184 for these Z values. For 
Z values much lower than those considered in this work, 
a spherical shell gap at A^ = 184 can not be ruled out 
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within this model. Comprehensive calculations for much 
of the superdeformed region with this model in future 
would clarify the situation. 

A lowering in energy of a large number of Nilsson or- 
bitals in going from ^^^115 to •^°''121 is accompanied by a 
slight increase in the quadrupole deformation. However, 
splitting of levels by the strong deformation leads to a 
significant increase in the binding energy of nuclei. A 
sudden jump in the /?2 value from the nucleus ^°°121 to 
•^'^'*123 also results in a kink in several levels and a further 
lowering in the energy of several Nilsson orbitals is antic- 
ipated. Therefore, this adds to the total binding energy 
of nuclei such as •^°'*123 and '^°^125. The consequence 
of this has already been seen in Figs. 5 and 6, where 
the Qa value of heavier superheavies is smaller by ^ 1-2 
MeV as compared to FRDM. This would inevitably pro- 
vide a larger stability and consequently a larger a-decay 
half-lives as compared to the FRDM. Thus, the present 
calculations with NL-SVl present an optimistic picture 
for producing superheavy nuclei heavier than Z — 115. 



E. Shell effects in nuclei and consequences 

As mentioned in Section I, the spin-orbit interaction 
and consequently its effect on the shell gaps and shell 
effects has a major influence on creation of magic num- 
bers. It is, however, still not clear as to how the spin- 
orbit interaction extrapolates in the extreme limits of the 
periodic table such as near the r-process path and for 
extreme superheavy nuclei. We believe that both these 
regions have a bearing upon each other. Inevitably, how 
the shell effects would transcend and extrapolate in un- 
known regions would be decisive in r-process nucleosyn- 
thesis of heavy nuclei as well as in synthesis of extreme 
superheavy nuclei. We put up the case in Section II that 
if the shell effects are strong along the line of stability, 
these would extrapolate strongly in the extreme regions 
and vice-versa. This was demonstrated reasonably well 
in ref. Q- 

Here a comment on the strength of the shell effects 
would be in order. It is well known that the macroscopic- 
microscopic approaches such as FRDM exhibit shell 
effects that are stronger at iV = 82 and N = 132 in go- 
ing to extremely rich regions of the r-process path. The 
strong shell effects at these magic numbers in FRDM 
has a consequence on the r-process nuclear abundances. 
Using the data from the FRDM in network chain calcu- 
lations, there is a considerable shortfall (troughs) in the 
abundance of r-process nucleosynthesis of nuclei below 
the peaks at A - 130 and A - 190 [111- These peaks 
in the abundance curve correspond to the above neu- 
tron magic numbers. Therefore, the strong shell effects 
in FRDM are not consistent with the r-process nuclear 
abundances which require a weakening of the shell effects 
in going to the r-process region. The strong nature of the 
shell effects in the FRDM is evidently carried away also 
in the region of superheavy nuclei. This becomes appar- 



ent in going to extreme superheavies in the neighborhood 
of Z = 120 as shown in Fig. 5 and 6, where large values 
of Qa are predicted by the FRDM. 

The macroscopic-microscopic calculations of ref. |23| 
are of the similar origin as the FRDM. These calculations 
have been successful in describing superheavy nuclei be- 
low Z = 109 as illustrated in Fig. 2 of ref. As 
mentioned earlier, results of ref. overestimate the 

Qa values (see Fig. 5 and 6) for superheavy nuclei above 
Z = 111 jlSj . However, comparing the general trend 
of the predictions from macroscopic-microscopic calcula- 
tions of ref. 01 with those of the FRDM near Z = 120, 
one finds a reasonably good agreement between the two 
model predictions on Qa values. Both the models predict 
larger values for Qa in the extreme region. This suggests 
that similar to the FRDM, strong shell effects are also 
carried over to the extreme region in the calculations of 
ref. 22]. This would indeed lead to very small a-decay 
half-lives for very heavy superheavy nuclei. On the other 
hand, the less strong shell effects in microscopic calcu- 
lations with NL-SVl predict Qa values that are smaller 
than those of the FRDM. This shows an immense im- 
portance of as to how the shell effects extrapolate in the 
extreme regions and how the shell effects influence prop- 
erties of nuclei in these regions. 



VI. CONCLUSION 

We have performed relativistic mean-field calculations 
for the recently observed alpha-decay chains of nuclei 
288115 and ^s^ns using the Lagrangian set NL-SVl with 
the vector self-coupling of w-meson. The Qa values 
have been calculated. It is shown that across the ex- 
perimentally observed chains, the NL-SVl values show a 
rather smooth behaviour with the change in the Z val- 
ues of the nuclei. The Qa values obtained with NL- 
SVl without blocking of odd-particle pairing describe 
much of the experimental data very well. In compari- 
son, the macroscopic-microscopic calculations |22j over- 
estimate the experimental values systematically. The a- 
decay half-lives calculated with Qa values from NL-SVl 
using the Viola-Seaborg systematics show a reasonably 
good agreement with the experimental values deduced 
from the recent experiment 11^ with the exception of a 
few nuclei. 

The decay products of the both the chains with Z — 
107, 109, 111 and 113 are shown to be prolate deformed 
with a deformation varying from (32 ~ 0.22 for Z = 107 
to P2 0.17 for Z = 113. Curiously, the nuclei ^^^m 
and 28^113 in both the chains exhibit a shape-coexistence 
of a low prolate shape with (32 ~ 0.17 and a superde- 
formed prolate shape with P2 ^ 0.50. In both the cases, 
the low (32 shape corresponds to the lowest energy min- 
imum. The parent nuclei ^^^115 and ^^^115 of both the 
chains, however, exhibit a very large deformation akin to 
a superdeformed shape. 

We have also made predictions for heavier elements in 
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the two chains of 115 and ^^'115. It is shown that with 
NL-SVl the ground-state of the higher Z counterparts of 
both the chains are superdeformed with /?2 values in the 
vicinity of 0.60 with a tendency of an increase in fii value 
in going to heavier elements. This can be attributed to 
a softening of the shell effects in going to the extreme 
case. The Qa values obtained with NL-SVl for nuclei 
with Z = 117, 119, 121, 123 and 125 are about 1-2 MeV 
lower than the predictions of the FRDM This has 
a consequence that a-decay half-lives of the heavier nu- 
clei are 2-3 orders of magnitude larger than the existing 
predictions of very short half-lives from the macroscopic- 



microscopic models |23, |5^. This scenario would put 
synthesis of heavier superheavy elements on a feasible 
footing. 
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